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Abstract 

Electrons obeying the Dirac equation are investigated under the non-relativistic c H > oo 
limit. General solutions are given by derivatives of relativistic invariant functions which 
possess discontinuity at the light-cone, yielding the delta function of (cf) 2 — x 2 . This 
light-cone singularity does survive in this limit to show that the charge and the current 
densities of electrons travel with the speed of light in spite of their massiveness. 


1 Introduction 

It is well known that the Dirac electron has a piece vibrating with the light velocity, called 
” Zitterbewegung” whose origin is supposed as a mixing of the positive and the negative 
components. Discussions have been made by the Heisenberg equation of motion for the 
Dirac Hamiltonian [TJ or by the momentum representation for the Dirac field 
Electron moving in the zigzag motion with the light speed also appears on the stage of 
the pilot-wave approach to quantum field theory [3] (Feynman had already discussed the 
zigzag motion with the light velocity in the context of path integral[4])- 

We shall, in this paper, focus the electron field \P(x) itself, to show that electrons 
obeying the Dirac equation inevitably bear portions traveling with the light speed in the 
non-relativistic limit. As a preliminary, let us recall the Dirac Hamiltonian 

Hd = c (a • p + f3mc ) ; a = 7°7 ; /3 = 7 0 , ( 1 ) 

with 

= diagGm = (1, —1, —1, —1) ; ([i,v = 0,1,2,3) , (2) 

being the 4x4 gamma matrices represented as 

<x: Pauli matrices . (3) 
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7°=( l _°! I >T = 


0 a 

— a 0 



In order to investigate the non-relativistic limit, it is useful to perform a unitary (called 
the Foldy-Wouthuysen) transformation [5], 


U = exp 


IP 

. me 


e(p) 




(4) 


with 


such that 


tan 2 —9(p) 
me 


\P\ 

me 


UHvU ] = pH ; H = c^fp^mP? 


(5) 


Here the lower component corresponds to the negative energy state which should be 
discarded in the non-relativistic world. Therefore we shall pick up the positive energy 
part and study its wave mechanical structure in the next section. In the following section 
[3J we shall treat covariant solutions of the free Dirac equation, which will be extended 
to interacting cases in sec|H The final section is devoted to the discussion. Some of the 
detailed calculations in secQ]are relegated to the appendices. 


2 Wave Mechanics of H = cy 7 p 2 + m 2 c 2 

Consider a single component wave mechanics governed by the Hamiltonian ([5]) , that is, a 
wave function obeying the Schrodinger equation 


The solution reads 


x ) = cy / yhVp _ Trn^c2 <F(f, x) . 




( 6 ) 


(7) 


where ?/>(0, x') is an arbitrary function, 


K(Ax\t) = (a:| exp 
f d 3 p 


—i—\/P 2 + m 2 c 2 
h 


\x 


(2irhy 


exp 


— (p ■ Ax — ct\/p 2 + m*c 


2„2 


Ax = x — x 


( 8 ) 


is the kerneQ with P designating the momentum operator and the integration range from 
—oo to oo has been omitted here and hereafter unless otherwise specified. 

In order to calculate the kernel (JHJ) , introduce k = p/h and write 

me 


P 


h 


(9) 


Hn the path integral formalism, this is a typical example that the Hamiltonian prescription is more general [B] 
than Feynman’s [7], whose Euclidean case, that is, t —it in J5]) has been discussed in [3]. 
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to obtain 


r ri 3 k 

K(Ax-t) = j ^3 exp 


ik ■ Ax — ictyjk 2 + jj? 


I(r,t) = / dkex p i (kr — cty/ k 2 + /r 2 J 

./ fVl ' 


1 dl(r , t) 
47 r 2 r dr 

r = |Ae| , 


( 10 ) 

( 11 ) 


where use has been made of the polar coordinates to the final expression. Put k = n sinh 0 
to find 

coo 


/ oo 

dO cosh 0 exp[— i/i (ct cosh 0 — r sinh 0 )] 

-OO 

whose exponent reads, with the aid of an addition theorem, as 


( 12 ) 


ct cosh 0 — r sinh 0 = < 


x 2 cosh (0 — a) : tanh a = — ; ct > r 


ct 

ct 


j - ££ 

— x 2 sinh(0 — /3 ): tanh /3 = — ; r > ct 


(13) 


with 


xl = (ct) 2 - r 5 


(14) 


Make sifts 0 i—/ 0 + a, /? and again utilize the addition theorem to obtain 

/ OO 

d0 cosh 0e~ lx+ cosh0 


+/r#(r — ct) 


• — OO 

coo 


/ oo /*oo 

d0 cosh 0e ix ~ sinh 0 + sinh p dB sinh 0e ix ~ sinh 0 

-OO J —OO 


with 


X± = C\ ±x 2 u 


(15) 


(16) 


/ OO 

d0 sinh 0e~* x+ cosh 0 = 0.) By noting 

-OO 

it 

/ OO /‘OO n ,- 

d0 cosh 0e ix ~ sinh0 = / c?(sinh 0 )e* x_sinh 0 = 2v rd( X -) = -SU-x 2 ) , 

-oo J —oo C * 


that 


and cosh/3 = r/,,/—a ; 2 in view of (|13l) . the second term of (|15li reads 




fi6(r — ct) cosh p I d0 cosh 0 e* x_ smh0 = 0(r — ct) I n 


2nr 


—xf 


2 7T 


C 


— 6 (\h x 


= Anr6(r — ct)5(x z fl ) = 6(r — ct) ——5(ci — r) = m5(ct — r) 
where we have used the relations, 

Kyf-x 2 ) = 2 \[-^A x V> > 6 '(°) = \ ' 


(17) 
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Finally the Bessel function formulas [9] 


/ oo 

d& cosh 0 e~ ix+ cosh e 

-OO 

/ OO 

dO sinh O e ix ~ sinh 0 

-OO 


-M 2 \x + ) 

2iKi(x~) , 


lead us to 


I(r,t) 


ir5(ct 


r ) — 2ctp 2 


9(ct — r) 


^7r fff 2) (x + ) ^ 


9(r — ct) [i 


Mx-) 


X- 


(18) 


Here note that I consists of different functions in the regions ct > r and ct < r, which 
causes the delta function singularity when a differentiation is made. (The //-independent 
delta function emerges from the huge momentum domain p me.) 

The kernel is, from () 10 D . obtained, by differentiating (1181) with respect to r, as 


K(Ax ; t) = 

ctp 2 
' 2n 2 


1 d 

47 rr dr 


5(ct — r) 


ip 

4-/r 2 


5(ct — r) 


9(ct-r )| 



i9(r 



where use has been made of 


d_ 

dr 


9{ct — r) 


—5(ct — r) ; 



ct) 


5{ct — r) , 


and then (see Appendix lAl) 


ctp 2 
2iT 2 r 


5{ct — r) 


■x H?\x+) , .. ^(x-) 

2 X+ X- 


-h S{f *~r) ’ 


as well as [TO] . 

zdz \ z J 



J n ,N n , (also H^ 2) 


) and K n , 


(19) 


( 20 ) 


( 21 ) 


by noting 


1 d HU 9 1 d 

-7T- = T/i 2 -W- 

r dr x± °X± 


( 22 ) 


By taking the non-relativistic limit c >->■ oo, that is, p e-)- oo ([9]), the kernel (fT9l) reads 

K(Ax\t) / '=°° ~ \^ x \) + o(^p 3/2 ^j , 

where we have employed the asymptotic expansion of the Bessel function j 111 


(23) 


Z n (z) z =°° OfV 1 / 2 ) , : J n ,N n , (also H&) ; K n (z) z =°° ^=(1 + 0(z~ 2 )) . (24) 
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Now put 


il>(0,x') = ( —\ 
an ) 


3/4 


/ 2 ' 


exp 


X 

2 a 


then substitute (1231) into ([7|) with changing the variables as x' H > Ax to find 


&(t,x) = - 


r 

/ d 3 (Ax)5 (ct — \Ax\) 


iy? / 1 ^ 3/4 
47T \an / 

which, by introducing the polar coordinates, yields 

(ct — |*|) 2 


exp 


(x + Ax)‘ 
2 a 


+ 0(/r 3 / 2 ) , (25) 


&(t,x) = - 


iy 2 a 1 ^ 4 ct 


IT 


3/4 


\X\ 


exp 


2a 


exp 


(ct + \x\f 
2 a 


When ct, \x\ S> y/a, the second term fades away and around the peak, 

ct — |a?| ~ 0 , 

becomes 

(ct — |aj|) 2 ' 


&(t,x) w - 




IT 


3/4 


■ exp 


2 a 


(26) 


(27) 


(28) 


which apparently travels with the speed of light in spite of its massiveness. 

The origin lies in the delta function 5(ct — |ai|) (1191) emerging from the discontinuity 
of I (1181) between the time-, ct > \x\, and the space-like, ct < |a:|, regions. We shall call 
this as a light-cone singularity. (It should be emphasized that the light-cone singularity 
cannot become visible under the momentum representation [13] . ) 


As a necessary consequence, any wave written as 

^(x) = I d 4 yV(x,y)i’(y) ; v V>(y) , 

must have the light-cone singularity, if T> contains derivatives to some function with a 
discontinuity on the light-cone. The solutions of the Dirac equation meet with this, so, in 
the next sections, we shall study those. 


3 The Charge and the Current Density of Free 
Electrons 

First let us summarize the relativistic invariant functions which participate in solving a 
relativistic equation. The D-dimensional scalar Klein-Gordon field is given as 

(i n + /J ?)<Hx) = j(x ), n = d»d„ = -A-'r,-A, (29) 

j=l 3 

where y is defined by (|9|) and J(x) is a source, complicated function of cj> describing 
interactions. Here and hereafter the repeated indices always imply the summation. When 
J = 0, the solution is 

4>o( x ) = j d D yA(x - y)ip(y) , (30) 
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where A(x) is an invariant function defined by 


(□ + y 2 )A(x) = 0 , 

A(x) = - J d D ke{k 0 )5{k 2 - y 2 )e~ lkx 

= ~ (2 J , (31) 

and ip(y) is an arbitrary function. The notations, 

k 2 = (ko) 2 — k 2 ; kx = koxo — k ■ x , (xq = ct) , 


with k and x being (D — l)-dimensional vector, should be understood. When J / 0, the 
solution is 

<t>{x) = 4>o{x) ~ J d D yA F (x - y)J{y ) , (32) 

where 4>o(x) is (f30l) and A F (x) is the Feynman propagator, 


(□ + y 2 )A F (x) = -S D (x ) (= -(5(x 0 )5(xi) • • • <5 (xd-i)) , 

e“ 


(33) 


A f (x) = J 


d D k "~ ikx 


(2n) D k 2 — y 2 + ie 


These are shown as m 


A{x)=- ^ )u O{x 2 )e{x 0 ) f^yVx 2 ^ J- v [yVx 


D- 2 

v = —tt— > (34) 


Af(x) = - 


k 


2v 


0(x 2 ) (^yVx 2 ^ | J-u^yyfx 2 ^ — iN_ u (^yVx 


4(2 ir) u 

+i 'j 6(—x 2 ) (ji'J—x 2 ^ K_ 


—ar 


(35) 


with x 2 = (x 0 ) 2 -x 2 . By noting J- n (z) = (-) n J n (z), N_ n (z) = ( -) n N n (z ) and K_ n (z) = 
K n (z ), they become in D = &{v = 1)[I5], 


A(x) = 
A f (x) 


e(xo) 

2 7T 


5(x 2 ) - f —0{x z ) 


/U 2 2 .Ji(yVx 2 ) 


1 xf 2 \ , k 2 0{x 2 ) J^yVx 2 ) 


—%- 


.y 2 0(—x 2 ) Ki(fj,V—x 2 ) 


4-7T 2 


yV^ 


X* 


. N\ (yVx 2 ) 


(36) 


(37) 


Note that they have a discontinuity on the light-cone. (Any relativistic invariant function 
does.) In view of ()30p and (1321) . however, there are no derivative so that we cannot have 
light-cone singularities for scalar fieldfU. (We do not care the 0(1) delta function in (1361) 
and (137(1 under the non-relativistic limit y t—)■ oo.) 

2 As for complex scalars, the current density J^x) = ie {(j)* {x)d l _ l (j){x) — (d^cj)* (x))<f>(x)) does have a light-cone 
singularity. 
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The free electron field &o(x) obeys the Dirac equation, 

(ifi - y)%{x) = 0 , ^ = ^9^, 

where y is @ and <Po(x) is given as 

<Po{x) \ 


Vn(x) = 


Xo(x) 


( 38 ) 


(39) 


with (po and xo being two component spinor. 

The solution of (l38l) reads, 

Vote) = J d A yS{x - y)ip(y) , 

where S(x) is the invariant function for the Dirac field, 

{if) — y)S{x) = 0 , S(x) = {if) + y)A{x) 


(40) 


(41) 


with A{x) (13611 and ijj{y) being an arbitrary four component spinor. From (1401) and (1411) . 
<Po(x) must own the light-cone singularity. 

For the sake of simplicity, the initial electron configuration is assumed to be, 


1 %) = S(yo)f(y) f q° 


where £o is a constant two component spinor. Since from (13111 (with D = 4) 
A(x) = 0 ; <9q A(x) = —5 3 {x) , 


xo=0 


x o =0 


(14011 with (I42H implies the initial condition, 

%{x 0 = 0,x) = —if(y) 
In the following we consider three cases; 


£o 

0 


(42) 


(43) 


(44) 


f(y) 


f {3) (y) = 


l 


rV4 


, v 

exp - — 


2 a 


/< 2 | M^% S )e*p(-g) ; 

f {l \y) = a 1/4 d{y 2 )5{y 3 )exp 


2 _ 2 , 2 
2/2 = 2 / 1 + 2/2 


(45) 


v\ 

2 a 


which is called as 3-, 2-, and 1-dinrensional packets respectivehU. 
In view of (l4H) . the solution (1401) becomes, 




<p ( 0 k \x) 

X ( 0 k) (x ) 


= d 3 y 


id* + ,i 

-icr • V x 


\A(x 0 ,x - y)f (k) {y)£, 0 ; (A: = 1,2,3). (46) 


3 In terms of the length scale L, the dimension of a spinor field ^{x) is L 3 / 2 . A(x) is L 2 then S(x) is L 3 
from (l4ll) . Hence i/j(x) (l40l) is L -5 / 2 so that f(x) l l42l) is L -3 / 2 . 
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In the following, we shall discuss the charge and the current densities defined bjQ 

(fc)_ e , Tr (k)\ Tr (k) ef ( fe )t (fc) (fc)t ( fc ) A 

pi = — iPq = ~ (¥>o ( x m (*)+ X 0 (*)Xo ( x )) ’ 

7 (t) = eS^* 1 Vi '= e Vo** V)* 7 *? V) + XoAxVvVV)) . (47) 


ZS(x) " 3 ” 


Now take the non-relativistic limit c e-)- oo (/j, i->- oo) to find that Zl(x) (f36l) reduces to 

(48) 


(49) 


4vr u,\i x 

whose derivatives are calculated as follows: first note 

d 0 (e(x o )0(x 2 )) = 2 \xq\5(x 2 ) , 

by use of 

<9 0 e(x 0 ) = 25(a:o) ; 5(xq)B(x 2 ) = 0 ; x 0 e(x 0 ) = |x 0 | 

as well as do = 2xod/dx 2 and 

d 


dx 2 


9(±x z ) = ±5(x 2 ) . 


Thus 


a A , y liV 2^ J l(I i ^ I ) , B 2 , \ c\l 2no d ( 

doA(x) = —2|x 0 | 5{x )- j== -h — e(xo)0(x )2x 0 7tft I - ^=— I 


47r 


whose second term reads 


47T 


dx 2 




2 d (JitflVx 2 )^ z= l iVx s 2 d ( 

^ dx 2 \ pVx 1 y ^ dx 2 dz\ z ) 


(50) 


® = o(n 3/2 ) , 


where J n (z ) ~ =°° 0(1/y/z) (1241) has been considered in the final expression. Meanwhile 
the first term becomes, with the aid of J\(z)/z\ Q = 1/2 (j!28j) . yielding to 

d 0 A(x) = l ^\xo\5(x 2 ) + o{^ 3,2 ^j . (51) 

Similarly by dk = —2 Xkd/dx 2 and (1501) 

d k A(x) =-j-x k e(x 0 )S(x 2 ) + o(^n 3/2 ^J . (52) 

From these, we can convince that the leading terms in the non-relativistic approximation 
are nothing but the light-cone singularities. Therefore by noting that 


5(x 2 ) = 


2\x\ 


< 5 (x 0 - |*|) + < 5 (*o + 1 * 1 ) 


(53) 


4 In relativity, a current 4-vector reads = ( cp , J) in the MKS unit. The Dirac particle \V induces a current 

jm = with W = ^y 0 . 















as well as pA(x) = 0(/i 3 / 2 ), (1461) reads as 


= 


l -^-jd 3 y5(x 0 - \y\) ^ -y ■ a j f {k) (y + *)Co + 0(V /2 ) , 


(54) 


where we have noticed xq > 0 and made a shift y i—>• y + x . 

Now proceed to individual cases: from (1451) and (1541) the 3-dimensional solution is 


8-7T a 3 / 4 


* 2 \ 

exp i — — I 


2 a J 


1 

-o ■ V 


-f (3) (*o,*)Co , 


where 


I {3 \x 0 ,x) = J d 3 y5(x 0 -\y \) 


x 0 


exp 


* • y jr_ 

a 2 a 


which becomes with the aid of the polar coordinates 


/(3) (x0);k) = ^£ 

\x\ 


exp 


2 a 


exp 


|*|*o 


exp 


\X\Xq 


Hence 


( 3 ) _ ip 2 a 1 / 4 x 0 


= 


4 \x\ 


1 - 


a exp 


*0 * 


1 


1 + 


a exp 


xo\x\ 

When xq, |*| 2> y/a, we have 


(*o + |*| f 

2 a 


(*o - |*|) 2 

2 a 

( Co 

V n-<r£o 


x 


n = 


* 


% 


( 3 ) _ ip 2 a 1 ^x 0 
4| * | 


exp 


which further turns out to be 

(3) ip 2 a 1 / 4 

*Pq ' ~---exp 


(*o - 1*1 f 

2 a 


(*o ~ |*|) 2 

2 a 


Co 

n ■ er£ 0 


Co 

n ■ o £o 


around the peak, 


*o — |*| ~ 0 , 


which apparently travels with the speed of light. 
Insert (f59l) into (f47]l to find 


(55) 


(56) 


(57) 


T~"T ' ( 58 ) 


(59) 


(60) 


p^ = 


2e ( / u 2 a 1 / 4 


exp 


(*o - |*|) 2 ^ f 


cko, j (3) = 


cnp 


(3) 


( 61 ) 
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where use has been made of the anti-commutation relations {aj, a*,} = 25jk in Those 
travel with the speed of light. 

It would be easier to prepare packets restricted in the 2- or 1-dimensional region. The 
former reads, by inserting (H5j) into (JMJ), as 


( 2 ) = ifl 2 (_*l_ 

0 87raV4 P 1 2a 


1 


-(72 


d / (2) (x 0 ,X2)£o • 


(62) 


x 0 dx 2 


where er 2 = (cri,cr 2 ),x 2 = (xi,x 2 ) and 

/ ( 2 ) (x 0 ,x 2 ) = j d 2 y 2 5(x 0 - \y 2 \) 


exp 


x 2 ~ V2 V2r_ 
2 a 


a 


(63) 


Here we have put X 3 = 0, since the observation should also be made in the X 3 = 0 plane. 
(1631) becomes under the polar coordinates, to 


n 00 r2n 

I (2 \x 0 ,x 2 )=/ dyy d65(x 0 -y) 
Jo Jo 

^ r, 

2 a/ 7 0 


exp 


l* 2 |y . y 2 

-cos #- 

a 2 a 


= xq exp 


d6 exp ( — 


\x 2 \x 0 


cos 6 


When | x 2 1, xq S> \/a, the saddle point 6 = ir gives a asymptotic value such that 


I { 2 \x 0 ,x 2 ) ps J 2 ^° exp ^ 


2 a a 


(64) 


Thus (1621) reads 


( 2 > VV/4 
0 P 


(xq ^ \X 2 \f 

2 a 


£0 

^2 • CT2^0 


(65) 


Here the final expressions has been obtained by putting xo/|* 2 1 1 in the coefficient; 

since the peak is now given as 

xq - \x 2 \ ss 0 , ( 66 ) 

whose velocity is again the light speed. The charge and the current densities (I47|) read 

^ - 7 ( • J(2) - MM ■ (-) 


Finally we discuss the 1-dimensional case written, after putting x 2 = 0, X 3 = 0, as 

1 


( 1 ) */r 2 a 1 / 4 

<Frj =-exp , 

0 8 vr 2 a 


a_ d | / ( 1 ) (x 0 ,xi )^ 0 , 

x 0 


( 68 ) 


where 


/ ( 1 ) (x 0 ,x 1 ) = / dyi5(x 0 - |?/i|) 


exp 


_ m 

a 2 a 


(69) 
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which turns out, by recalling xq > 0, to be 


7 (1) (x 0 ,xi) = 


dyi5(x 0 + yi)+ dyi5(x 0 -yi) 
> Jo 


exp 


x\y\ _ yf 

a 2 a 


exp 


( XqXi 
\ a 


)+-(~) 


exp 


2a 


Thus 


^(i) = iy 2 a 1/4 
0 8vr 


exp 


(x 0 -xi) 2 \ / (x 0 + xi ) 2 

1 ± exp - 


2a 


When xq,x\ -^/a, it becomes 


(i) iy?a 1 / 4 

^ ' « ———exp 


8-7T 


2 a 


(xp - xi) 2 \ / £o 
2a / \cri£o 


£o 

- V°i£o. 


whose peak is 

x 0 - xi « 0 , 

around which the charge and the current densities (1471) read 

,2 


P (1) = 


2e /V« 1/4N \ f (x 0 -a:i) 2 ^ t fl) / a) n n \ 

) eXP (-a- ) ^ ’ J ’°’°) 


(70) 

(71) 

(72) 

(73) 

(74) 


In view of (1611) . (1671) and (I74|) . the free electrons travel with the speed of light in the 
non-relativistic world!. 


4 Electrons in a Laboratory 

In a realistic situation, electrons interact with the electromagnetic field A^(x) such that 

(ip ~ y)^(x) = -^~'y IJ 'Af J ,(x)W r (x) , (75) 

whose solution is 

&(x) = &o(x) + J d 4 yS F (x - y)e^A /x (y)^(y) , (76) 

where &o( x ) is the free field discussed in the previous section and S F (x) is the Feynman 
propagator for the Dirac held, 


(iP ~ v)Sf(x) = 6 4 (x) , S F (x) = (ip + h)A f (x) , (77) 

with 4p(s) given in ([37]) . 

The interaction is assumed to take place for a finite interval, giving 

= ft (fc, (y) ( o ) . (* = 1 . 2 , 3 ), (78) 

5 Note that the relation between the current and the charge density J = cnp in (ED- (ED, and m is merely 
a kinematical consequence by taking the upper component £0 only. 
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where 3-, 2-, and 1-dimensional packets are 

1 


h { 3 ) (y 0 ,y) = 
^ 2 \yo,y) = 
h { 1 \y 0 ,y) = 


(r + 6) 1 / 4 r 1 / 4 6 3 / 4 

1 

(r + 6) 1 / 4 r 1 / 4 6 1 / 4 
6 1 / 4 


exp 


ysL^yl 

' 2t 2b 


vl vl 


% 3 )exp(-g-g); yl^yl+yl, (79) 


(r + 6) 1 / 4 r 1 / 4 


^( 2 / 2 )<5( 2 / 3 ) exp ( -|^ - || ) , 


with £1 being a constant two component spinor. (Note that the dimension of h(y) is 
L~ 5 / 2 . See the footnote [3j) If we write 




<Ph ] ( x ) 

Xh\x) 


= / d A y 


id 0 + y 

-itr ■ V x 


|^f(*-2/)/i (&) (2/)6 ; (* = 1,2,3), (80) 


the general solution (1761) reads as 

^(fc) = + 5 

with given by (176[) . 

In the non-relativistic limit y i->- 00 , 4p(s) (l37[) reduces to 


(81) 


Z\ F (^ 


/iM -00 y 2 6{x 2 ) 
“ 8ir 


Ji (yy/x 2 ) . Ni (yy/x 2 ) 


— 1- 


— *- 


.y 2 0(—x 2 ) K\(yy/ —x 2 ) 


4vr 2 


ar 


y 2 6(x 2 ) H[ 2 ' (yVx 2 ) .y?6{—x 2 ) K\ (yy/—x 2 ) 


8vr n\/x^ uy/ -i‘ 


47T 2 


(82) 


where we have utilized (11251) to the final expression whose form reminds us (1181) in sec 12J 
The derivative reads 


d fl A F (x) = 2x fl —5(x 2 ) 

o 7 T 


y 2 2\ JiiyVx 2 ) 2i (it Ni (yVx 2 ) Ki(yy/—x 2 ) 


7 r \ 2 


yV—3 


+ o(/// 2 ) 


where use has been made of ([501) and the asymptotic behavior (f24l) . Noting (j!28 j) and 
(11311) in Appendix |A] we have the light-cone singularity, 


d^A F {x) = x^^dix 2 ) + 0(V /2 ) 


(83) 


Contrary to the previous situation, we now need both terms in ([53]) . since yo cannot 
always be positive, so that 

2 

(ido + y)A F (x -y) = S(x 0 -yo~\x- y\)-5{x 0 -y 0 + \x- y\) + 0(V /2 ) , (84) 

■ __.■, , . -y 2 cr-(x — y) 

- ia -V x A F (x-y)=-^ - 

87r^|x — y | 

Then (1501) turns out to be 


S(x 0 -y 0 -\x - y\)+5(x 0 -y 0 +\x - y\) +o[y 3/2S j . (85) 


< k) - wfr 


x 0 


-\x-y\,y)Th( k \x 0 + \x-y\,y)') -y • a Ki , (86) 
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where a shift y i—)■ y + * has been made. 

Let us exam individual cases: with the 3-dimensional packet 

1 y /4 

h 8ir 2 \ t 6 3 (t + 6 )/ 


exp 


x\ 


0 


* 


2 r 2b 


rd o 

bcr ■ V 


read, 


Ih\ x o,*)£i , 


where 


4 3) (x 0 ,*) = 


= 2 


d 3 y | 

d 3 y 


exp 


xq\v\ 


+ exp - 


xo\y\ 


exp 


/ 


\y\ 


■ exp 


(r + b)y 2 
2rb 


cosh 


3=0 , 


2/1 exp 


(r + b)y 2 
2rb 
x • y 


x • y 


which, with the aid of the polar coordinates and the error function) 16]. 

2 r x 

erf (a:) = —p= / dfexp(—f 2 ) , 

V n Jo 

yields to 

(i>*o + t |*|) 2 


4 3) (*o,*) = 


(27r6) 3 r 1 
t + b |a:| 

—erf 


erf 


bx o + r|*| 


\J2 rb(r + b) 

( bx o — r |*| \ 
v / 2r6(r + 6) 


exp 


exp 


2 r6(r + 6) 

(6*o - r|*|) 2 
2 t6(t + 6) 


Applying the differentiations, in (1571) and assuming *o> 1*1 ^ Vd we obtain 


V, 


(3) 


—y 2 ( rb 3 




2\/2tt \(t + 6) 3 


1/4 


exp 


(x 0 - |*[f 

2 (r + 6) 


6 

n ■ cr^i 


around the peak (1601) . (Details are relegated to Appendix iBl) 

Therefore in view of #^ (1551) an( l 4.^ fl2H) • the total ^(jUJ) is obtained as 


„• ..2 „l/4 
^(3) = a 


exp 


r + a + b 


(*o - |*|) : 


where we have introduced a two component spinor 

t — a + b, , 2 


^ = exp - 


4a(r + b) 
rb 3 


4 a(r + 6) 

>onent s 

(*o - |*|)^ ) £o 


£ 

n • er£ 


^ w 7T \a(r + b ) 3 
The charge and the current densities are 


1/4 


exp 


f t — a + b 
\4a(r + b) 


(* 0 - l*l) 2 6 


( 3 ) = £^( 3 )V( 3 ) = f UlEl 


2 «.v 4 


Ptot 


exp 


t + a + b 
2a(r + 6) 


(*o - l*D 2 ^) 


2e( y 2 a 1 / A 


,2 r 


exp 




a 

r6 3 


t6 3 


7r V a(r + b ) 3 


exp 


(*o - |*|) 2 ^ t 


2(r + b) 


(87) 


( 88 ) 


(89) 


(90) 


(91) 


(92) 


(93) 


£6 


'2 / rb 3 \ 7 ( r + a + b 

+ nl n{wTW) exp {~M^Tb) 


(*0- I*|) 2 ) (cki -£ko) 


J t (3) ee e^ 3 ) f 7 V (3) = cnpg . 


(94) 
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Since each term has a peak at the light-cone, the electron signal traveling with the light 
speed would be observed. 

Next we consider the 2-dimensional case: again we restrict x to be X 3 = 0 so that (18611 
with (179ll are found as 


af’ = 


-/* 


87t 2 \t6(t + 6) 

where 

Ih ] ( x 0,252) = 2 J 

which yields 


1/4 


d 2 y 2 

12/1 


exp 


exp 


2 2 \ / 

' r<9| 

*0 * 2 ) 


2r 2b ) \ 

&<x 2 • - 

(x + b) 2 

*2 • 1/2 

0,7. ^2 



d_ 

dx 2 


Ih\ x 0 ’X 2)6 , 


cosh 




4 2) (*o,a!2) 


(27r) 2 r& 2 


\x 2 \{t\x 2 \ + 6x 0 ) 

+ 


exp 


(27r) 2 r6 2 


*21(1-1x21 - 6X 0 ) 


(r|x 2 | + fr*o) S 

2 r6(r + 6) 

'( r|x 2 | - 6x 0 ) 2 
2 rb(r + 6) 


exp 


Therefore when xo, |® 2 | v^x- \/& 

1/4 




( 2 ) 


-y 


rb 3 


x|x 2 | + bx 0 


exp 


(*o - l*2| 


47t V(r + 6)5 J V |x 2 | 1 2(r + b) 

(Details are relegated to Appendix [Bj) This further reduces to 


6 

n 2 • <x 2 £ 1 




( 2 ) 


-V 


47t \(t + b) 3 


T 6 s y /4 

exp 


(X 0 - |X 2 |) 2 

2 (r + 6) 


around the peak (l66ll . 

From (1651) and (l99ll . the total jF^ dHTTl is given by 


I ,( 2 ) = H^ a ^ f (r + a + b) 


Ay/2 7T 

e w =expf- 


(xo - |x 2 | 


\ 4a(r + b ) 

(x 0 - |x 2 |) 2 J^0 

(x — a + b) 


( 2 )__/ (x-a + 6),_ n 2 


4a(r + 6) 

2 / rb 3 xV4 


Cl N 

■ 0- 2 £i 


' ^(2) > 
n 2 • <t 2 £ (2) 


(x 0 - |x 2 |) 2 6 


7r v a ( r + b) 3 J ^ V 4a(r + 6) 

The charge and the current densities are thus found as 

& - *( wH - 


2e/ /i 2 a 1 / 4 V 


(95) 


(96) 


(97) 


(98) 


(99) 


( 100 ) 


( 101 ) 


- 


y 4:V2^ J 

+i 


exp 


(x 0 - |x 2 |) 2 ^ ctc , 2 

?ns0 + 


7r \a(r + b) 3 


2 a 

tP \ 1/4 

exp 


r5 3 


7r \a(r + 6) 3 

(r + a + b) 

2 a(r + 6) 


1/2 

exp 


(X 0 - |X 2 |) 2 ^ , 


2(r -I- 6) 

(x 0 - |*2|) 2 ") -Cko) 




t( 2 ) — 
•'tot — 


= e!F (2)t 7 0 7 lF (2) = (cn 2 p { S t , 0 


( 102 ) 
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Finally we consider the 1-dimensional case: by restricting x to be X 2 = x% = 0 ([79]) 
brings (186)) to 


^ 


8tt 2 \r(r + b) 
dyi 
\yi\ 


1/4 

exp 


x o x i 
2b 


J\b(Tidi J 


Ih\x 0,®l)£l 


I^\x 0 ,xi) = J 


exp 


(r + 6) 2 

2 Tb Vl 


xm\ 

b J 


2 cosh 




After a little calculation (see Appendix [B]) it reads, when xq,x\ ^ Vr, Vb, 


V, 


(i) 


— p? ( rh 3 


1/4 

exp 


(x 0 - X! f 


h 4 ^ 2 T\(r + by 
From (HID and (11051) the total ^ 1 - 1 (IHTT) is obtained as 

(r + a + b) 


2(r + b ) J l cri£i 


6 


,■,.21/4 
^(1) = ^ 


87T 


exp 


4a(r + b) 


■(xo-XiY 


ew S exp(- ^;.7_ u ^ (xo-x,) 2 )^ 


(r - a + 6) 


4a(r + 6) 
+ n/I 


r6 3 


7t \a(r + ft) 3 
The charge and the current densities arcjfl 


1/4 


exp 


() 


(r — a + b) 
4 a(r + b ) 


(x 0 -a:i) 2 6 


(b = £^(i)t^(i) = 2e ^ 


2 nV 4 


Ptot 


exp 


2e 

C i 87T 


exp 


87t I 
(xq-xi) 2 


(r + o + 6) 
2 a(r + 6) 


(XO-S!) 2 C 


di) t Ai) 


7T 




7T \a(r + 6) 3 


2a 

rt 3 \ 1/4 

exp 


Tb* V /2 
exp 


^ 0+ 2 Va(r + 6) 3 
(r + a + b) 


(103) 

(104) 


(105) 


(106) 


(107) 


(xp-xi) 2 ^ 


2 a(r + 6) 


2(r + 6) 

(*o - ^i) 2 ) (&i - eko) 


e{6 


7"! 1 ) 

'tot 


= e^ 1 ) t 7 °7 = (cpS,0,0) 


(108) 


5 Discussion 

In this paper, first we discuss the wave mechanics of H = c\Jp 2 + m 2 c 2 , which tells us 
that solutions of the Schrodinger equation inevitably possess the light-speed portion in 
the non-relativistic limit c i —> oo. The reason is that the kernel contains a derivative acting 
to a function which owns the discontinuity on the light-cone. The solutions of the Dirac 
equation are also expressed by differentiations to the invariant functions A(x) and A?{x) 
which consist of different functions in the time- and the space-like region, thus yield the 
light-cone singularity, which was the contents of seci3]and[4l In relativistic field theories, 

6 Again note that J = cnp in (l94l) . (11021) . and (11081) does not imply the current’s velocity is c. See the 
footnote [5] 
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c appears as p = mc/fr, so that the non-relativistic limit implies /j t-> oo which, however, 
also interprets the semiclassical h i —> 0 or an infinite mass limit According to 

the last case we can convince ourselves of survival of the light-cone singularity for massive 
particles. We should emphasize that our conclusion has been derived exclusively in the 
^-representation of <F(x) not in the momentum representation. 

The situation is unchanged if the source (1781) would have a velocity v: consider, for 
example, 


h {i \y) ^ d 3 (y - Py 0 )exp {; (3 = ^ 


Then from 


l h. ~ jd 4 y5(x o - y 0 =F \x - y\)S 3 (y - /3y 0 )exp^~ 
= Jdy 0 5(x 0 -y 0 ^\x- /3y 0 |)exp^- 


yj 

2 r 

2 

0 

2 r 


(109) 


( 110 ) 


(The spinor part is irrelevant.) Since the zeros in the delta function are given by 


xo ~ 2/o = ±|* ~ Pyo\ 


y o = 


(*o + |®|) (1 - P ■ n) + 0(/3 2 ) 
(xo - |*|) (1 + P-n) + 0(/3 2 ) 


it reads when xq , * 2> \fr 


d'h 



(xq - |*|) 2 (1 + p ■ n) 2 \ 
2r ) 


( 111 ) 


which again shows that the maximum signal travels with the speed of light. 

In order to widen the possibility we finally consider the case that the initial configu¬ 
ration is given with a definite momenta p, that is, instead of the packets (|42|) take 

,2 s 


V’p(y) = 


in the solution 


1 


(a7r) 3 / 4 


<%o)exp( -p • y - 


y_ 

2 a 


Co 

0 


Ide o = 1 


d'pix) = j d 4 yS(x - y)ip p (y) . 

Since at xo = 0, as was in (144|) . (I112D implies an initial configuration, 

„ 2 \ 


dpi'Xo = 0, x) = 


— X 


exp 


% ar 

— p ■ X - 

h 2a J 


(an) 3 / 4 

whose momentum reads as 

d 3 x ^p(x) (— iKV ) d' p (x) = p 
Following a similar procedure from (1461) to (156 j) . we find 


Co 

0 


/■ 


d'p(x) = 


ip 


87r(a7r) 3 / 4 \h 


exp —p ■ x — 


Xq + x 2 



dp Co 


( 112 ) 


(113) 


(114) 


(115) 


(116) 
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with 


I P = d 3 y5(x 0 -\y\)exp[^-p-y- 


x ■ y 


which is further rewritten as 

I p = exp/ d 3 y5(x 0 - |y|) 

ia 


= 27raxo exp I — —p • V -—7 


\x\ 


with the help of (1571) . By noting 


V = n 


d 

d\x\ 


exp 


exp 


n = 


x • y 

a 

\x\x 0 

a 


x 


\x\ 


exp 


\x\x 0 


and that exp( — ~yP ■ n , ) is a shift operator, it reads 


h d\x 
exp 


Ip — ItTCLXQ' 


x\xq p ■ n 
~~a 1 h~ X ° I exp l 


\x\x 0 .pn 

-t-z—^*o 

a h 


\x\ — iap ■ n/h 


Therefore 


zpV /4 x 0 

p ' 47t 3 / 4 \x\ — iap ■ n/h 


7 


1 


\ xq(\x\ — iap ■ n/h) 

{ -1 


exp 


- 




i(p • nxo — p ■ x) (xq — \x\) 2 

2 a 


h 


+ 


1 + 


exp 


/ i(p ■ nxQ - p ■ x) (xq + \x\) 2 \ 
\ h 2a ) 


\ xo(|#| — iap ■ n/h) t 

When xq, |*| S> yfa ; \x\ — iap ■ n/h ss |*|, it yields to a plane wave, 


(117) 


(118) 


(119) 


( 120 ) 


Co 

ya • n£ 0 , 


.( 121 ) 


<£„(x) 


zp 2 a 4 / 4 
47r 3 / 4 


■ exp 


i(p ■ nx 0 — p ■ x) 


h 


Co 

t cr • n£ 0y 


( 122 ) 


around the peak xq — \x\ ~ 0 , which implies that the energy-momentum relation is given 

by 


E = c\p\ . (123) 

Therefore an alternative way to observe the relativistic remnant is a measurement for the 
energy and the momentum of electrons in the vacuum. 
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A The Derivation of ( 1201) 

Owing to the delta function S(ct — r) 

X± = vV ± (( ci ) 2 - r 2 ) ^ 0 , 


so that 


vr h [ 2 \ X+ ) 


X+ 


x +=0 X- 


7T Jl(X+) 


X-=0 

_ . fvriVi(x + ) 


X+=0 


2 X+ 


^i(X-) 


x+=o X- 


2 X+ 
where we have noticed 

h£\z) = J n (z) - iN n (z ) , (n = 0,1,2,...) 
From the definition of J n and I n m 


X - 0 


■«*> = (!)"£ 
k=0 

«*> = (?)”£ 


<-)' 


2fc 


k=0 


we have 


While 


•Mx+) 


X+ 


x+ 


=o 2 


^ Ni(x+) 


2 X+ 
tfi(X-) 


X- 


= Ji(x+) 
x+=o x+ 

= h(X-) 
x -=0 X- 


fc!(n + k)\ \2J 

5 

1 

fZ\ 2k 

k\(n + k)l \2) 

5 

h(x-) 


1 

X- 

X-=0 

2 ’ 

n X+ ) 

1 

1 

2 ) 

x+=o 4 

x+ 

X \ 

1 

1 


H - 7 

H- o - 

2 / 

x-=o 4 

xi 


X+=0 

x-=° 


with Euler’s constant 7 , in view of nn 
; N n (z) = J n (z ) (7 + In 


7T 

2 _ 


2k 


I HE 

'2 \ 2 / ^ fc!(n + /c)! V 2 

1 / 2 \ n y^ (n — k — 1 )! (z\ 2 fc 

"2I1J ^ at 

v 7 fc =0 


n+/c 


E 1 + E 1 

z ' m z —' m 

_ra=l m=l 


K»(z) = (-)" + 1 « 2 ) (7 + I 17 ) 


(-)” ^ 


5 £ 


(-)* 


k=0 


k\(n + A:)! \2 


2k 


k n+k 

E i + E 1 

Z —' m z —' m 

,m= 1 ra=l 


+ 


1 / 2 \ n ^H (n — k — 1 )! 


2 \ z 


E 

k=0 


k\ 


(124) 

(125) 

(126) 

(127) 

(128) 


(129) 


(130) 
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then 


k Ni(x+) 
2 X+ 

1 


Ki(x~) 


x+=o 


X- 


X •- 0 


= - (In x x 2 — In A / — 




x2=° 


H 2 xl /i 2 x 2 


x2=° 


4 L 


ln(ct — r) + ln(ct + r) — ln(r — ct ) — ln(r + ct ) 


ct=r 




where we have used ln(r — ct) = in + ln(ci — r) to the final expression. 
Therefore (11241) reads 


* H[ 2) (x+) 

2 X+ 


x+=o x- 


7T 7T l l 

x—=o _ 4 _ 4 + 2 _ 2 


□ 


(132) 


B Derivations of and 


(1901) : With the aid of the polar coordinates, (1881) becomes 

(r + % 2 


4 3) (*o , x ) = [ dy exp 

1*1 Jo 


2 rb 


cosh 


xo \ . , (\x\ 

TV smh (T 


iirbT, ( (r + % 2 

— / dy exp - 

x Jo \ 


2rb 


sinh 


bx o + r|*| 
rb 


y ) — sinh 


bx o — t\x\ 


rb 


-y 


,(133) 


where the addition theorem for the hyperbolic function has been used. In view of the 
error function formula (1891) we have 



sinh By 


1 

2 





so that (11331) becomes (TOD . 


(19T1) > : In view of (l90l) and (|87l) note 


exp 


x%_x?_ (top ±r|a;|) 2 \ = / (x 0 t|*1) 2 ^ 

2r 2b 2Tb(r + b) ) \ 2(r + b) ) 


(134) 


(135) 
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then apply (% and V to obtain 
(3), x_ -M 2 ( i~b 3 V /4 

^ ( ' 3 ) 2V2^V(^ + fe) 7 J 


bx o + r|®| 


erf 


® 


/ frxp + t|®| \ / (x 0 

\ yj2rb(r + b) J V 2 ( 


- ® 


2 (r + 6) 


6xq — r|®| 


erf 


® 


bxQ — t\ x\ 
^2rb{r + b ) 


| exp 


(x 0 + |®[) ; 

2(r + 6) 


X/i u 2V2^u^ + ^) 7 y 


1/4 


I'bxo + T \x\ b{ T + b)\^( bxo + T\x 
\ \x\ X 2 


+ 


+ 


bx o — t\x\ b{r + b) 


\x\ 


X A 


^ erf f 


Iv^ + b ), 

bx 0 — r \x\ 

^ y / 2r6(r + b) 


exp 


exp 


(X 0 - \x\f 

2 (r + 6) 

( x 0 + |®|) 2 
2(t + 6) 


2 rb(r + 6) 2 


7 r cc 


In ® 

eXP ‘-27"26 


(n • <r)Ci 


When xq, \x\ ^/r, V^, terms vanish, except the first one in the right hand side to give 




(3) 


/ rb 3 \^bxo + r\x 


h ~ 2V2^\(r + by) 

f bx o + 


xerf 


y a/2t6(t + b) 


\x\ 


exp 


which further, by noting 


(xp - \x\f 
2 (t + b) 

bx o + t|® 


& 

n ■ aii 


lirn erf(X) = 1 ; X= , _ 

-X^i ^-oo ^J2rb(r + b) 


becomes 


( 3 ) -y 2 ( rb 3 V /4 bx o + x \ x \ 

h ~2 V2^W + b)\ 


\x\ 


exp 


Oo - \x\y 

2 (r + b) 


6 

n ■ aii 


yielding to (19T1) around the peak (1601) . 

(j97|) : (l96jl reads by use of the polar coordinate 


r(2) 


(x 0 , ® 2 ) =2, J dy exp ^ V 2 ^ cosh ( J dO exp 


\ x 2\y 


cos 6 


When | X 2 | Vb, the saddle point method around 9 = n brings the angular part to 


2nb f\x 2 \y 
\x 2 \y 6XP V b 


1 + 0 


l*2| 


Then 


I h\ x o,x 2 ) 


2irb 


\x 2 \ Jo Vv 


exp 


(r + b ) y 2 + (r|® 2 | + bx p) 


2 rb 


rb 


, , (r + b) 2 (r|® 2 1 - bx 0 ) 

+ exp — 7n—y h- 7 - y 


2 rb 


rb 
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which is further rewritten, in terms of a dimensionless quantity Y, 


t\x 2 \ ± bx 0 
V = —^ Y 


as 


I^\ x 0 ,* 2 ) = 


_p(±) = 


't|ce 2 | ± &xq Z 100 dY 


f 


Vrt Jo VY eXP \ (rb) 2 [2 r + 6' y/ 

Since |x 2 |,xo S> \/r, the saddle point method around y = \/rb/{r + 6) gives 

{t\x 2 \ ± bx 0 ) 2 ' 


(r + b)(r\x 2 \ ± 6xq) 2 f y 2 \/r& 


( 136 ) 


y 


_p(±) = 


27 Trb 


■ exp 


t\x 2 \ ± bx o V 2t6(t + 6) 

Inserting this into (113611 we have (|97j) . 


1 + 0 


(t|* 2 | ± &x 0 ) 2 


( < Y\ 

(1551) w, : Apply differentiations in (l95l) and note the relation (11351) by putting x i->- x 2 


to obtain 


'A 


( 2 ) 




r6 3 


47t \(r + 6) 5 


1/4 


(r|* 2 | + 6x 0 


t\ x 2 \ - bx 0 
+ W-i—i-exp 

V *2 


where we have omitted terms of 

r6(r + b ) 


1*21 

( x 0 + |* 2 |) 2 

2(r + b ) 


■ exp 


(Xp - 1X2^ 

2 (r + 5) 

a ^ 

n 2 • 0- 2 £i ) 


(137) 


O 


O 


(t|x 2 | ± bxo) 2 J 

Under xo, 1*21 ^ \/t, \/&, (I137|) becomes 
(11051) d f h l> : (I104D becomes (we have put y\ y) 


b(r + b) 


|* 2 |(t|* 2 | ± bx o) 


!+„,*.) = 2 j T|exp(--h±V 


) cosh ( 


6xq + rxi 


y + cosh 


f bx o 


— TXl 


-y 


.(138) 


2 rb ) V rb ) \ rb 

Apply the differentiation in (|103l) , with the aid of the error function formula (1134j) , to find 


<A (1) (x 0 ,xi) = 


-y 


rb 3 


1/4 


, ( bx 0 + TXl 

erf Ut(r + h) 1 exp 


4y/2ir V( r + b) 3 

, c ( bx o - txi \ ( 

(2rtyT5)) exp ( 

which yields, when xq, xi 2> yjr, Vb, to (11051) . 


(xp + x{y 
2 (r + b) 


(x 0 - Xif 

2 (r + 6) 

f a > 
v <ua / 


( 139 ) 
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